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We  report  on  the  effects  of  different  multifilamentary  geometries  on  the  dynamic  current  and  held 
distributions  in  YBa2Cu307_^  thin  films.  Finite-element  model  simulations  are  presented  that  allow 
a  direct  comparison  of  the  calculated  magnetic  flux  and  current  profiles  with  experimental  data 
obtained  by  time-resolved  magneto- optical  imaging.  We  find  that  increasing  the  distance  between 
filaments  reduces  the  magnetic  coupling  and  decreases  the  ac  losses.  Screening  and  transport 
currents  also  redistribute  more  evenly  between  ah  filaments.  The  measured  transport  current  density 
profiles  are  in  good  agreement  with  the  modified  critical  state  model.  ©  2009  American  Institute  of 
Physics,  [doi:  10. 1063/1. 3212976] 


I.  INTRODUCTION 

The  separation  of  high-temperature  superconducting 
(HTS)  tapes  into  filaments  is  a  viable  approach  to  reduce  ac 
and  hysteretic  losses  in  high-power  applications,  where  alter¬ 
nating  currents  and  magnetic  fields  are  present 
simultaneously.1  Methods  such  as  mechanical2  or  laser 
scribing,3  photolithography,4  or  direct  printing  on  buffered 
substrates  using  inkjet  deposition5  have  been  used  to  create 
the  filamentary  structure  in  second-generation  HTS  coated 
conductors.  However,  ac  losses  of  finely  striated  tapes  can 
still  be  larger  than  desired  as  predicted  by  analytical 
expressions.6  Even  so,  the  magnetic  held  concentration  in  the 
nonsuperconducting  trenches  influences  the  flux  distribution 
inside  the  superconducting  filaments,  and  this  magnetic  cou¬ 
pling  ultimately  creates  losses  that  depend  on  the  distance 
between  filaments.  This  is  a  complex  nonlinear  problem  that 
necessitates  an  in-depth  understanding  of  the  flux  and  current 
dynamics  in  multifilamentary  HTS  structures  as  both  depend 
strongly  on  temperature  and  history  of  the  applied  held  and 
transport  current. 

In  this  paper  we  report  on  the  effects  of  different  mul- 
tihlamentary  geometries  on  the  dynamic  current  and  held 
distributions  in  YBa2Cu307_(5  (YBCO)  thin  films.  We  present 
finite-element  model  (FEM)  simulations  that  allow  a  direct 
comparison  of  the  calculated  magnetic  flux  and  current  pro¬ 
files  with  experimental  data  obtained  by  time-resolved 
magneto-optical  imaging  (TRMOI).  The  comparison  pro¬ 
vides  both  a  test  of  the  model  itself  as  well  as  new  informa¬ 
tion  on  the  complex  behavior  of  YBCO  thin  films  observed 
previously.7  We  find  for  filamentary  samples  where  the  ends 
of  the  filaments  on  each  side  are  connected  via  the  supercon¬ 
ducting  him,  see  Fig.  1,  that  increasing  the  spacing  reduces 
the  magnetic  coupling  between  hlaments,  which  respond 
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more  independently  to  an  applied  held  and  ac.  Screening  and 
transport  currents  also  redistribute  more  evenly  between  ah 
hlaments.  The  reduced  magnetic  coupling  affects  held  and 
current  dynamics  and  decreases  ac  losses. 


II.  MODELING 

In  the  FEM  simulation  we  neglect  thermal  aspects  and 
focus  on  the  electromagnetic  behavior.  The  superconducting 
him  is  considered  to  be  at  a  constant  operating  temperature 
below  the  critical  temperature,  and  no  quenching  occurs.  In 
addition,  since  the  samples  are  long  and  straight,  we  utilize  a 
two-dimensional  model  considering  only  the  superconduct¬ 
ing  cross  section.  The  electromagnetic  behavior  is  described 
by  Faraday’s  law 


The  superconducting  property  is  incorporated  by  means  of  a 
nonlinear  relation  between  electric  held  and  current  density 


FIG.  1.  (Color  online)  Example  of  a  filamentary  YBCO  thin  film  (blue  area) 
grown  on  an  insulating  substrate  (pink)  with  metal  contact  pads  (green).  The 
dc  field  B  is  applied  perpendicularly  to  the  film.  The  screening  current  Jsc 
flows  circularly  in  the  sample,  while  the  ac  transport  current  Jtr  flows  along 
the  filaments. 
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This  relation  describes  quite  well  the  experimental  current- 
voltage  characteristics  of  superconducting  samples.  The 
power  index  n  represents  the  steepness  of  the  transition  from 
superconducting  to  normal  state  and  is  indicative  of  the 
physical  phenomena  of  flux  flow  and  flux  creep.  In  high- 
temperature  and  type-II  superconductors,  magnetic  flux  lines 
enter  as  quantized  vortices  at  very  small  values  of  the  mag¬ 
netic  field  Bcl.  For  this  reason  in  the  simulation  we  express 
the  B-H  relation  by  using  B  =  jm0H.9 

The  model  is  implemented  in  COMSOL’s  general  partial 
differential  equation  (PDE)  module  and  uses  two  magnetic 
field  components  as  state  variables.  First-order  edge  ele¬ 
ments  are  used  to  satisfy  the  zero-divergence  equation  for  the 
magnetic  field.  Details  about  the  model  implementation  are 
given  in  Ref.  10.  The  external  magnetic  field  is  incorporated 
by  using  Dirichlet  boundary  conditions  in  the  air  domain. 
The  transport  current  is  included  by  using  integral  con¬ 
straints  in  COMSOL’s  point  settings.  The  integral  constraints 
introduce  coupling  between  filaments  corresponding  to  the 
superconductive  region  connecting  the  filaments  in  the 
sample  as  depicted  in  Fig.  1. 

In  the  FEM  simulations,  we  consider  a  system  consisting 
of  the  air  domain  and  six  superconducting  filaments  con¬ 
nected  at  the  ends.  The  dimensions  chosen  for  the  system 
correspond  to  YBCO  thin  films  grown  by  pulsed  laser  depo¬ 
sition  on  a  LaA103  or  SrTi03  substrate  and  patterned  by 
photolithography  to  obtain  a  multifilamentary  structure  (Fig. 
1).  Since  the  substrate  has  very  poor  conductivity  and  is 
nonmagnetic,  it  has  been  neglected  in  the  simulated  geom¬ 
etry.  These  films  have  been  previously  studied  by  TRMOI 
under  the  same  experimental  conditions  used  in  the 
simulations.11 

The  model  is  not  restricted  to  any  number  of  filaments, 
dimensions,  or  aspect  ratios.  However,  a  high  aspect  ratio  of 
the  thin  film  geometry  introduces  a  large  number  of  nodes 
that  can  severely  affect  the  computation  time  especially  for 
systems  with  a  large  number  of  filaments.  The  simulation  of 
two  cycles  takes  about  3-4  h  on  a  workstation  equipped  with 
a  2.4  GHz  processor  and  3  Gbytes  of  random  access 
memory.  In  order  to  avoid  transient  effects,  three  complete 
cycles  of  the  current  in  50  steps  are  simulated.  An  example 
of  the  simulated  geometry  and  corresponding  mesh  is  de¬ 
picted  in  Fig.  2(a).  The  circular  outer  boundary  for  the  air 
domain  is  used  to  impose  Dirichlet  boundary  conditions.  The 
density  of  nodes  in  the  mesh  increases  near  the  filaments  as 
shown  in  the  detailed  view  of  Fig.  2(b). 

In  order  to  allow  direct  comparison  between  FEM  simu¬ 
lations  and  TRMOI  measurements  of  the  field  dynamics,  we 
calculate  the  magnetic  profiles  at  a  height  of  1  jum  over  the 
sample  surface  [see  Fig.  2(b)].  This  approach  has  two  impli¬ 
cations.  First,  it  reduces  slightly  the  intensity  of  the  simu¬ 
lated  field  profiles,  similarly  to  what  occurs  in  the  actual 
measurements,  and  second,  it  roughens  the  profiles  adding 
noise  due  to  the  coarsening  of  the  mesh.  The  latter  effect 
could  be  avoided  by  refining  the  mesh  in  that  region  but  at 
the  cost  of  an  increased  calculation  time. 


FIG.  2.  (Color  online)  (a)  Complete  mesh  and  (b)  a  detailed  view  near  the 
edge  of  a  superconducting  filament.  The  red  dashed  line  indicates  the  region 
of  the  simulated  held  profiles. 


III.  EXPERIMENTAL 

In  our  TRMOI  experiments,  the  sample  size  is  normally 
5X12  mm2  with  six  filaments.  The  striation  is  made  by  pho¬ 
tolithography  and  wet  etching.  The  average  width,  spacing, 
and  length  of  the  filaments  are  130  /zrr l,  200  jum,  and  4 
mm,  respectively.  Figure  1  displays  a  schematic  three- 
dimensional  drawing  of  the  samples  where  the  blue  areas 
indicate  the  superconducting  YBCO  film,  which  is  about  300 
nm  thick  grown  on  a  SrTi03  substrate.  The  samples  exhibit  a 
critical  transition  temperature  ( Tc )  of  91  K  with  a  sharp  tran¬ 
sition  within  2°  (at  2.2  Oe  loss  data)  as  determined  by  ac 
susceptibility  measurements.  The  samples  are  carefully  in¬ 
spected  for  defects,  which  may  couple  two  neighboring  fila¬ 
ments  (other  than  the  end  caps)  or  discontinue  the  current 
flow  in  a  filament. 

In  the  TRMOI  experiment,  the  sample  is  zero-field 
cooled  to  15  K  and  then  a  dc  magnetic  field  of  Ba  =  5  mT  is 
applied  perpendicular  to  the  YBCO  film  together  with  a  sinu¬ 
soidal  ac  transport  current  /0=  8  A  at  a  frequency  / 
=  1000  Hz.  The  TRMOI  experiment  is  described  in  detail  in 
Ref.  11.  Figure  3  displays  images  taken  from  a  six-filament 
YBCO  sample  at  different  phase  points  of  the  ac  transport 
current.  Bright  areas  in  the  image  correspond  to  regions  of 
high  magnetic  flux  density  near  the  edges  of  the  filaments. 
The  dark  areas  in  the  center  of  the  image  indicate  flux-free 
regions.  The  uniformly  gray  region  (red  in  the  online  ver¬ 
sion)  represents  the  intermediate  magnetic  field  outside  the 
sample.  Positive  and  symmetric  flux  penetrations  occur  on 
both  sides  of  the  sample  due  to  the  applied  dc  field  (Fig.  3, 
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FIG.  3.  (Color  online)  TRMOI  images  from  a  six-filament  YBCO  thin  film 
for  different  phase  points  of  the  applied  ac  with  70=8  A.  A  dc  magnetic  held 
of  5  mT  is  applied  perpendicular  to  the  him.  The  arrow  indicates  the  current 
direction. 


middle  image).  Additional  flux  penetrates  at  the  top  edge  of 
the  filaments  when  a  transport  current  flows  in  the  positive 
direction  (Fig.  3,  left  image),  while  the  same  effect  occurs  at 
the  bottom  edge  when  the  current  flows  in  the  negative  di¬ 
rection  (Fig.  3,  right  image).  This  phenomenon  is  most  pro¬ 
nounced  for  the  outer  filaments. 

IV.  RESULTS 

Figure  4  (top)  shows  the  magnetic  flux  distribution  cali¬ 
brated  from  light  intensity  profiles  of  the  TRMOI  measure- 


FIG.  4.  (Color  online)  The  magnetic  held  distribution  (top),  the  total  current 
distribution  (middle),  and  the  transport  current  distribution  (bottom)  of  the 
six-hlament  YBCO  sample  with  applied  ac  (70=8  A)  and  perpendicular  dc 
magnetic  held  (Ba= 5  mT).  The  positive-to-negative  ac  cycles  are  colored  in 
red-black-blue  shades.  The  vertical  dotted  lines  indicate  the  edges  of  each 
hlament. 


ments.  The  magnetic  flux  changes  its  sign  and  amplitude 
with  the  phase  evolution  of  the  ac.  The  positive  flux  has  a 
higher  absolute  value  than  the  negative  flux  because  the  ex¬ 
ternal  field  is  static  and  always  contributes  a  positive  flux. 
The  total  current  distributions  are  shown  in  Fig.  4  (middle) 
obtained  by  numerical  inversion  of  Biot-Savart  law  from  the 
magnetic  field  distributions.  The  total  current  flowing  in  the 
first  three  filaments  1-3  is  mostly  positive  while  it  is  pre¬ 
dominantly  negative  in  filaments  4-6.  For  all  51  phase- 
dependent  profiles,  the  current  distribution  outside  of  the  six 
superconducting  filaments  is  set  to  zero,  which  has  only  a 
small  effect  on  the  total  current  density  as  shown  further 
below.  The  transport  current  distributions  shown  in  Fig.  4 
(bottom)  are  calculated  by  subtracting  the  shielding  current 
determined  from  the  current  distribution  at  zero  phase  point. 
As  expected,  in  the  positive  (negative)  half-period,  only 
positive  (negative)  current  (red  shades)  flows  through  all  six 
filaments. 

Figure  5  shows  the  total  transport  current  flowing  in  all 
six  filaments  at  each  phase  point  (red  dots).  The  current  os¬ 
cillates  between  —7.72  and  8.33  A  (black  dashed  line), 
which  is  in  good  agreement  with  the  experimental  value  70 
=  8  A.  This  result  shows  that  the  current  distributions  ob¬ 
tained  from  the  TRMOI  measurements  are  quite  accurate. 
For  the  multifilamentary  film,  when  the  filaments  are  not 
coupled,  the  current  will  evenly  distribute  between  all  fila¬ 
ments.  On  the  other  hand,  when  the  filaments  are  coupled, 
the  current  flows  across  the  transverse  resistance  between 
filaments12  for  resistive  interfilamentary  connections  or,  in 
this  case,  via  superconducting  connections  at  the  ends.  Al¬ 
though  these  two  different  types  of  connection  (supercon¬ 
ducting  and  nonsuperconducting)  will  have  a  different  effect 
on  the  current  distribution,  they  both  will  result  in  nonuni¬ 
form  current  distribution  across  the  filaments.  Figure  6  dis¬ 
plays  the  transport  current  distribution  (solid  lines)  at  two 
extreme  phase  points,  </>=  90°  and  270°,  which  clearly  re¬ 
veals  that  the  six  filaments  are  coupled  in  our  experiment. 


V.  COMPARISON  OF  MODEL  TO  EXPERIMENT 

To  further  analyze  the  current  distribution,  we  use  the 
modified  critical  state  model  (CSM)  for  the  superconducting 
thin  film  geometry  with  an  applied  ac  and  a  static  magnetic 
field  applied  perpendicular  to  the  sample  face.11  The  blue  and 
red  dashed  lines  are  fits  of  the  modified  CSM  for  a  single- 
bridged  sample,  which  indicates  the  overall  current  density 
distribution  across  all  filaments.  The  green  dotted  curves  rep¬ 
resent  fits  of  individual  filaments,  which  agree  fairly  well 
with  the  measured  data  in  shape  and  amplitude.  Compared  to 
the  transport  current  values  calculated  from  the  measured 
data,  the  average  error  is  approximately  15%. 

We  also  compare  the  TRMOI  data  with  FEM  simulation 
data.  Figure  7  shows  the  magnetic  field  distribution  (top)  and 
the  total  current  distribution  (bottom)  obtained  from  FEM 
calculations.  The  parameters  for  the  FEM  simulation  are  the 
same  as  in  the  measurement  Ba- 5  mT  and  70=8  A  at  1000 
Hz.  However,  the  critical  current  density,  /c=4.1 
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FIG.  5.  (Color  online)  The  integrated  transport  current  (red  dots)  in  the  six-filament  YBCO  sample  and  best  fit  to  a  sinusoidal  curve  (dashed  black  line). 


X  107  A  cm-2,  is  more  than  twice  the  value  of  the  single- 
bridged  YBCO  thin  film.  The  increase  in  Jc  is  unexpected 
and  requires  further  investigation. 

Figure  8  shows  a  comparison  between  FEM  simulations 
(dashed  lines)  and  the  TRMOI  data  (solid  lines)  for  the  mag¬ 
netic  field  distributions  (top)  and  the  total  current  density 
profiles  (bottom)  at  three  phase  points:  (fi=  90°  (blue),  (fi 
=  270°  (red),  and  (fi=  180°  (green).  The  profiles  match  quite 
well  for  the  field  distributions.  The  peaks  are  higher  in  the 
FEM  simulation  than  in  the  TRMOI  data.  Moreover,  the 
measured  field  penetrates  further  inside  of  the  superconduct¬ 
ing  filaments.  This  behavior  is  also  apparent  in  the  current 
density  profiles.  For  example,  in  the  center  of  the  first  fila¬ 
ment,  the  saddle  shape  total  current  density  profiles  have 
lower  values  for  the  FEM  simulations  than  the  measured 
data.  This  implies  that  the  Jc  used  in  the  FEM  simulation  is  a 


little  higher  than  in  the  experimental  data.  The  FEM  simula¬ 
tion  also  has  higher  absolute  peak  amplitudes  than  the  TR¬ 
MOI  data,  which  also  means  that  the  FEM  simulation  has 
slightly  higher  Jc.  Nevertheless,  the  data  from  the  multifila¬ 
mentary  sample  yield  a  critical  current  density  at  least  twice 
as  large  as  the  Jc  of  single-bridged  YBCO  thin  films,  which 
is  a  significant  advantage  of  the  multifilamentary  film. 

VI.  EFFECT  OF  FILAMENT  SPACING 

Next,  we  investigate  the  effect  of  magnetic  coupling  be¬ 
tween  superconducting  filaments  by  studying  the  dynamic 
current  and  field  distributions  as  a  function  of  distance  be¬ 
tween  filaments.  The  width  of  the  filaments  is  fixed  at 
120  jmm  while  the  spacing  is  varied.  We  simulated  three  dif¬ 
ferent  geometries  with  width/distance  ratios  8/1,  4/1,  and  2/1, 


FIG.  6.  (Color  online)  Comparison  between  measured  transport  current  density  profiles  (solid  lines)  and  fits  of  the  modified  CSM  model  (dotted  lines)  for  two 
phase  points:  90°  and  270°. 
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FIG.  7.  (Color  online)  FEM  simulation  of  the  magnetic  field  distribution 
(top)  and  the  total  current  distribution  (bottom)  for  the  six-filament  YBCO 
sample  with  applied  ac  (/0=8  A)  and  perpendicular  dc  magnetic  field  ( Ba 
=  5  mT).  The  positive-to-negative  ac  cycles  are  colored  in  red-yellow -blue 
shades.  The  vertical  dotted  lines  indicate  the  outer  edges  of  the  six  filaments. 

comparable  to  multifilamentary  coated  conductors  reported 
in  literature.2-6  The  applied  field  intensity  is  5  mT  and  the  ac 
I(t)=I0  sin(27 rft)  has  an  amplitude  70= 8  A  and  a  frequency 
f- 1000  Hz.  For  the  superconducting  material,  we  choose 
Jc= 3X1011  A/m2,  i^lO-4  V/m,  and  n- 25,  which  are 
typical  values  for  YBCO  thin  films  at  about  30  K. 

A.  Simulation 

In  Fig.  9  we  present  field  profiles  at  different  phase 
points  for  three  different  geometries  used  in  the  simulations. 
We  also  plot  magnetic  flux  distributions  as  time-dependent 
field  maps  in  Figs.  9(d)-9(f)  to  present  more  clearly  the  time 
variations  of  the  field.  The  magnetic  flux  partially  penetrates 
into  the  filaments  from  the  edges  forming  flux  fronts.  These 
regions,  as  well  as  the  flux-free  regions  (Meissner  state),  re¬ 
main  constant  during  the  cycle  in  all  the  filaments  for  the 
three  different  geometries  shown  in  Figs.  9(a)-9(c).  The  flux 
profiles  at  the  edges  of  the  filaments  alter  during  the  cycle 
due  to  the  self-field  of  the  applied  ac.  The  magnetic  flux 
distribution  changes  also  for  the  three  different  geometries. 
Increasing  the  filament  distance  allows  more  flux  to  pen¬ 
etrate  into  the  space  between  the  filaments.  Here,  the  peak 
intensity  of  the  field  profiles  grows  linearly  with  distance 
[Figs.  9(b)  and  9(c)].  The  flux  profiles  reach  their  maximum 


FIG.  8.  (Color  online)  Comparison  between  FEM  simulations  (dashed  lines) 
and  measured  magnetic  field  distribution  (top)  and  total  current  distribution 
(bottom)  for  three  phase  points:  </>= 90°  (blue),  <fi=  180°  (green),  and  </> 
=  270°  (red). 

value  of  65  mT  at  the  edges  of  the  two  outer  filaments.  The 
change  in  geometry  does  not  affect  the  maximum  value  of 
the  peaks  in  the  flux  profiles,  but  it  alters  the  dips  that  form 
at  the  external  edges  of  the  outer  filaments,  as  indicated  by 
black  arrows  in  Figs.  9(a)-9(c).  This  reveals  how  the  sample 
screens  the  external  magnetic  flux  for  the  different  geom¬ 
etries.  As  the  filament  spacing  increases  the  coupling  be¬ 
tween  filaments  is  reduced  and  the  screening  behavior  be¬ 
comes  more  similar  to  independent  filaments. 

The  current  dynamics  provides  additional  useful  infor¬ 
mation  regarding  the  coupling  between  filaments.  Figure  10 
shows  the  current  density  profiles  at  different  phase  points 
and  the  corresponding  time-dependent  field  maps.  The  maps 
show  at  first  glance  how  the  current  in  each  filament  changes 
sign  during  each  cycle,  passing  from  positive  (red)  through 
zero  (white)  to  negative  (blue).  Interestingly,  a  “return  cur¬ 
rent”  of  opposite  sign  [black  arrows  in  Figs.  10(a)-10(c)] 
flows  in  a  narrow  region  at  the  edges  of  the  left  outer  fila¬ 
ment  in  the  first  half  cycle  (-180°  <  <f><  0°)  and  on  the  right 
outer  filament  in  the  second  half  cycle  (0°  <  cf)<  180°).  This 
current  causes  the  steep  variation  from  positive  to  negative 
flux  observed  in  the  same  region  of  the  field  profiles  (Fig.  9). 
A  similar  behavior  has  been  observed  in  TRMOI  measure¬ 
ments  of  single-bridged  YBCO  thin  films.13  The  return  cur¬ 
rent  is  reduced  by  increasing  the  spacing  between  the  fila- 


Downloaded  17  Mar  2010  to  134.131.125.50.  Redistribution  subject  to  AIP  license  or  copyright;  see  http://jap.aip.org/jap/copyright.jsp 


063904-6  Lucarelli  et  al. 


J.  Appl.  Phys.  106,  063904  (2009) 


rmnrn 

EEE 


FIG.  9.  (Color  online)  Magnetic  field  profiles  at  different  phase  points  and  corresponding  time-dependent  intensity  maps  for  three  multifilamentary  geometries 
with  filament/groove  ratios  corresponding  to  8/1  [(a)  and  (d)],  4/1  [(b)  and  (e)],  and  2/1  [(c)  and  (f)]. 


ments  from  15  to  60  /xm.  The  increased  distance  reduces  the 
coupling  between  the  filaments  and  redistributes  the  screen¬ 
ing  currents  between  the  filaments.  The  growth  of  the  peak 
intensity  at  the  edges  of  all  the  filaments  reveals  the  reduced 
coupling  between  filaments.  For  the  three  different  geom¬ 
etries  the  current  density  in  the  two  outer  filaments  shows  the 
most  intense  value  corresponding  to  the  critical  current  den¬ 
sity  Jc=3  X  1011  A/m2,  as  predicted  by  the  CSM.  This  result 
yields  further  proof  of  the  accuracy  of  our  method. 


B.  Comparison  to  experiment 

The  current  density  profiles  shown  in  Fig.  10  are  ob¬ 
tained  from  the  field  profiles  by  inversion  of  Biot-Savart  law. 
This  method  has  some  limitations  but  it  allows  direct  com¬ 
parison  of  the  simulated  results  with  experimental  measure¬ 
ments.  The  current  profiles  obtained  by  the  inversion  method 
differ  only  slightly  from  those  calculated  directly  by  the 
COMSOL  model.  This  procedure  does  not  change  the  overall 
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FIG.  10.  (Color  online)  Current  density  profiles  at  different  phase  points  and  corresponding  time-dependent  intensity  maps  for  three  multifilamentary 
geometries  with  filament/groove  ratios  corresponding  to  8/1  [(a)  and  (d)],  4/1  [(b)  and  (e)],  and  2/1  [(c)  and  (f)]. 
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FIG.  11.  (Color  online)  Total  current  flowing  in  all  filaments  calculated  at 
different  phase  points  for  three  different  geometries  (symbols)  and  best  fit 
(solid  line). 

shape  of  the  current  profiles,  nor  does  it  affect  their  time- 
dependent  behavior,  but  it  smoothes  and  partially  reduces  the 
intensity  of  the  peaks. 

Figure  11  shows  the  integrated  current  densities  at  dif¬ 
ferent  phase  points  and  best  fits  for  the  three  geometries  used 
in  the  simulation.  The  screening  component  of  the  current  is 
symmetric,  and  thus  it  cancels  out  in  the  integration  of  the 
current  density  over  all  filaments,  and  we  obtain  the  transport 
current.  Figure  1 1  shows  that  the  amplitude  and  phase  of  the 
transport  current  are  not  affected  by  the  geometry.  This  quan¬ 
titative  analysis  also  confirms  that  the  reduction  in  current 
intensity  introduced  by  the  inversion  method  is  less  than 
1.3%. 

C.  ac  loss  consideration 

Figure  12  shows  the  transport  current  per  filament  for 
three  different  geometries.  Only  the  first  three  filaments  are 
represented  for  clarity.  The  behavior  of  the  other  three  fila¬ 
ments  is  completely  symmetric.  The  current  flows  in  each 
filament  without  any  significant  shift  in  phase  with  respect  to 
the  applied  current.  The  current  flows  mostly  in  the  outer 
filaments  for  all  geometries.  The  120/60  geometry  shows 
higher  values  of  the  current  in  the  internal  filaments  as  com¬ 
pared  to  the  other  geometries.  As  the  spacing  between  fila¬ 
ments  is  increased  the  transport  current  distributes  more 
evenly. 

We  also  use  the  FEM  to  calculate  the  electric  field.  From 
the  current  and  electric  field  profiles,  the  instantaneous 
power  losses  can  be  calculated  using  the  formula 


FIG.  12.  (Color  online)  Transport  current  per  filament  in  the  first  three 
filaments  (inset)  for  three  different  geometries. 


♦(deg) 

-180  -135  -90  -45  0  45  90  135  180 


time  (10*3sec) 

FIG.  13.  (Color  online)  Instantaneous  ac  losses  for  four  different 
geometries. 

P(t)  =  \j(t)-E(t)  \.  (3) 

Figure  13  shows  the  time-dependent  losses  for  four  different 
geometries  with  the  filament  width-distance:  120-15,  120-30, 
120-60,  and  120-90.  The  instantaneous  losses  decrease  with 
increasing  filament  spacing  in  the  range  <fi=- 36 °  to  </>= 90°, 
while  in  other  regions  the  losses  slightly  increase.  However, 
the  values  of  the  losses  integrated  over  the  whole  cycle  show 
an  effective  reduction  of  about  30%  in  the  120-90  geometry 
as  compared  to  120-15. 

VII.  CONCLUSIONS 

In  summary,  we  studied  dynamic  field  and  current  dis¬ 
tributions  in  multifilamentary  YBCO  thin  films  by  TRMOI 
and  by  numerical  simulation  using  COMSOL’s  PDE  module 
(general  form)  to  solve  Maxwell’s  equations.  A  comparison 
of  the  TRMOI  data  with  FEM  simulation  shows  that  the 
multifilamentary  films  have  at  least  twice  the  critical  current 
density  Jc  than  single-bridged  YBCO  thin  films.  Further¬ 
more,  an  increase  in  filament  spacing  alters  the  field  and 
current  dynamics  due  to  reduced  magnetic  coupling  between 
filaments.  Screening  and  transport  currents  redistribute  more 
evenly  between  all  filaments,  which  respond  more  indepen¬ 
dently.  It  also  significantly  reduces  the  losses  of  multifila¬ 
mentary  YBCO  thin  films. 
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